Degenerate ground state and quantum tunneling in rotating condensates 
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Quantum tunneling introduces a fundamental difference between classical and quantum mechanics. When- 
ever the classical ground state is non-unique (degenerate), quantum mechanics restore uniqueness thanks to 
tunneling. A condensate in a rotating trap with a vortex can have such a degenerate classical ground state, a 
degeneracy that is excluded in the absence of rotation at least when the Gross-Pitaevskii equation applies. If 
the rotating trap has a center of symmetry, like a figure eight (a peanut), the vortex may be on either side with 
the same energy yielding a degenerate ground state, a degeneracy lifted by quantum tunneling. We explain how 
to compute the rate of tunneling in the WKB limit by estimating the action of the trajectory in the Euclidean 
version of the dynamics. 



INTRODUCTION 

The search for physical "mesoscopic" situations where 
quantum tunneling could be significant is motivated by the 
search of quantum systems where off-diagonal correlations 
typical of quantum mechanics could be created and yield in- 
terferences without classical counterpart. We introduce be- 
low a novel physical system where such a tunneling exists. 
As well-known |1| Bose-Einstein condensates (BEC) rotat- 
ing fast enough have vortices in their ground state. We find 
that such vortices may have two different classical equilib- 
rium positions mapped into each other by geometrical sym- 
metry. Later we introduce the basic equations for describing 
quantum tunneling between the two degenerate states. This 
system is interesting because the Euclidean action allowing to 
estimate the tunneling frequency is complex, not purely imag- 
inary contrary to cases discussed in textbooks, and the single 
vortex of the "classical" state becomes, under the barrier, a 
pair of interacting vortices, something unnoticed before to the 
best of our knowledge. 

The wave-function common to all atoms in a BEC trapped 
in a potential F(r) is a solution of the Gross-Pitaveskii (GP) 
equation. This is a function in the ordinary sense because the 
quantum fluctuations can be neglected in the dilute limit, as 
shown long ago by Bogoliubov. Including in the G-P equation 
the possibility of constant rotation of vector one has: 

ih—- = /\^-ihYX^l-V^^g\^H^V{Y)^. (1) 

The trapping potential V(r) grows sufficiently fast at infin- 
ity, g is the coupling constant assumed positive and propor- 
tional to the s-wave scattering length, and m is the mass of 
each particle. A numerical version of equation ([T]) is found by 
taking as length scale ^ — , with n being the mean den- 
sity of particles over the surface a of the 2D condensate, 
n = ^ / |^(r)p(ir, and e a dimensionless parameter, arbi- 
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trary for the moment. The unit for ft is and the time 

unit is ^2^. With this system of units, equation (|TJ) becomes 

= -^A^ + iM^ + (I^P^ + V(r)^) , (2) 



where M is the real antisymmetric linear operator M = 
^ ~^^)' ^^^^ i^^y) being rectangular coordinates 
normal to the rotation axis and uj = \ft\ — ^ . Furthermore 
the potential V{t) has been made dimensionless. The ground 
state is a solution of ([T]) in the form of e~*^^^(r) with the 
smallest /i, for the given J dr|^p. Without rotation (ft = 0), 
the ground state is non-degenerate 1 1 1, but no such general re- 
sult exists with rotation, and we give below examples of the 
opposite. 

Consider a potential V(r) such that the condensate is con- 
strained to stay inside a closed curve shaped like a peanut. 
The idea behind this choice is that the vortex created by ro- 
tation may have an equilibrium position on either side of the 
peanut and so yield a degenerate ground state. We get such 
a curve from the Cassini oval, the locus of a point at which 
the product, 6^, of the distances to two foci of coordinates 
r = (±a, 0), is constant. The oval is made of two pieces for 
b < a and of a single one for 6 > a |2|. Physically the length 
a (for instance) is such that the dimensionless parameter e is 
small. The Cartesian equation of the oval is V{t) = with 

V(r) = + 4a^x^ - {x^ + + a^)^ 

Besides e, two dimensionless parameters characterize the sys- 
tem, one geometrical, the ratio ^ (assumed to be bigger than 
1), and the other measuring the strength of the perturbation 
due to rotation. Physically, this strength is found by compar- 
ing the circulation around a vortex and the one due to a solid 
body rotation at angular velocity uj. This yields the dimen- 
sionless ratio UJ = ^^9^^ , and defines e = — 7^=. One 
expects a transition for a finite value of uj from a vortex-free 
ground- state to a ground- state with vortex, a transition de- 
pending on the value of ^ . To find the ground state, we may 
simply consider the minimizers of the following nondimen- 
sionalized GP-functional, or energy: 

j {|V^p-2^(z^A^^) + ^|^|4 + ^|^p}dr, (3) 

where D is the domain defined by F(r) > and 5R(-) denotes 
the real part of its argument, ^ being the complex conjugate 



of ^. This minimum is constrained by the normaUzation con- 
dition 



/ \^\^dr= [ F+(r)dr, 
Jd J d 



where (r) denotes the positive part of the Cassini oval po- 
tential V^(r). Different solution branches are found for a range 
of the parameter values of a, 6, e and the rotation speed uj us- 
ing numerical methods similar to that described in 13J. We 
hereby present some examples. 

For a = 0.25, b = 0.275, and e = 0.001, Fig [T] displays the 
energy computed using ([3| for the various solution branches 
(i.e. the solutions without vortex, with one symmetric vortex, 
with one off-center vortex, with two symmetric vortices, and 
with three vortices) as a function of the rotation speed. The 
relative energy differences (1 — ^) are highlighted for uj near 
360, where En is the energy of a solution with at least one vor- 
tex, and £^0 the energy of the vortex free solution correspond- 
ing to the same rotation speed. Although the differences are 
very small, it is clear that in this range of uj, the solution with a 
non- symmetric (off-center) vortex has the lowest energy and 
makes the ground state. Therefore our initial guess is cor- 
rect: in this range of parameters values, the ground state with 
an off-center vortex is degenerate. The different solutions are 
shown in Fig|2] 
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FIG. 1: Energy values (left) and relative energy differences (right), 
as functions of the rotation speed, for different solutions with the 
Cassini oval potential. Black-square: no vortex. Red-circle: non- 
symmetric vortex. Blue- triangle: two vortices. Purple-diamond: 
three vortices, Green-star: one symmetric vortex. 

Similar energy diagrams can be found for different values 
of the parameters. With a smaller b = 0.27 but the same 
values of a and e, only three different solutions branches were 
found near = 415. We expect that a degenerate ground- state 
with a non- symmetric vortex exists in other geometries than 
the Cassini oval. We have replaced the Cassini oval potential 
F(r) by a piecewise constant potential of value +1 inside two 
squares connected through a rectangular channel, and (—10) 
elsewhere. When the channel width is that of the squares, 
D is 3. rectangle, and the ground state solutions keep the do- 
main symmetry. A non-symmetric vortex branch appears for 
a smaller channel width. We found that for a narrower chan- 
nel the energy differences of the different solution branches 





FIG. 2: Surface plots of |^| corresponding to the ground state with- 
out rotation (first row left, top view) and the five solutions found at 
the rotation speed uj = 360 (upside-down view, they are respectively 
solutions without any vortices in the interior oval, with a single non- 
symmetric vortex in one half of the oval, with a vortex in the center 
of the oval, with two symmetric vortices, and with three vortices). At 
this rotation speed, the ground state is the one with a non- symmetric 
vortex which is a degenerate state together with its mirror image with 
respect to the y-axis. 



are larger than for the Cassini oval when the ground state is 
degenerate (see Figs [3] and |4]). Even larger energy differences 
between the non- symmetric ground state and the other solu- 
tions are found by decreasing the potential in the channel con- 
necting the two sides. 
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FIG. 3: Energy values (left) and relative energy differences (right), as 
functions of the rotation speed, for different solutions with a piece- 
wise constant potential. Black-square: no vortex. Red-circle: non- 
symmetric vortex. Blue-triangle: two vortices. Purple-diamond: one 
symmetric vortex. 




FIG. 4: Plot of — 2/) I for a non-symmetric degenerate ground 
state with the piecewise constant potential at the rotation speed u = 
175. 

Based on our simulations, the following observations are in 
order: by increasing rotation speed with a standard symmetric 
harmonic potential, the ground state is given first LU by the 
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vortex free solution, then the single centrally symmetric vor- 
tex solution, then the solution with two and/or more vortices. 
It can be argued that, if in the range of the rotation speed where 
the ground state has only one vortex, there are substantial en- 
ergy differences between this ground state and other solutions 
(such as the vortex-free and the multi- vortex solutions), by 
properly perturbing such potentials in a non- symmetric fash- 
ion, non-symmetric ground state can be created and some can 
be degenerate. For instance, by lowering the potential away 
from the center (but not too far), symmetrically with respect to 
the ?/-axis, a couple of local pinning sites can be formed which 
would attract the central vortex to the vicinity of one of these 
sites, and lead to a solution whose energy is lower than the 
central vortex state. Yet if the perturbation is small, then the 
changes to the energies of the vortex free or multi- vortex solu- 
tions are also small, hence leaving the non-symmetric vortex 
as the new ground state. Such a solution can have the vortex 
situated at either one of the two pinning sites, thus becomes 
degenerate. This scenerio has been confirmed in our numeri- 
cal experiments as well. Other configurations of pinning sites, 
with three or four-fold symmetry for instance, could produce 
solutions with even higher degeneracy. The keys to realize this 
experimentally, aside for being able to produce the right geo- 
metric pattern, are to increase enough the pinning strength to 
differentiate the energy between the central vortex state and 
the pinned vortex state and to, at the same time, maintain 
substantial energy differences between the single vortex solu- 
tion and the other solution branches. Our simulations provide 
strong evidence to such possibilities. It did not escape our 
attention that a superconducting thin film in a perpendicular 
magnetic field yields a closely related situation, mathemati- 
cally speaking at least |4|. Therefore we expect that in certain 
range of parameters too the same kind of peanut like geometry 
yields there a degenerate ground state with a non-symmetric 
vortex. 

Whenever the classical ground-state is degenerate, one ex- 
pects quantum tunneling to restore the symmetry and make 
one single quantum ground- state out of the set of different 
classical ground-states derived from each other by symmetry. 
In the semi-classical approximation, the tunneling time Tq is 
proportional to exp(S'£;/fi), where 5'^; is the Euclidean clas- 
sical action, supposed to be much larger than h. If the system 
is very large, Tq may be so large that quantum tunneling is 
never observed. Nevertheless one may imagine intermediate 
cases, outside of the realm of atomic or nuclear phenomena, 
where the Euclidean action is not that large compared to h , 
leading to a non-negligible tunneling probability. The tunnel- 
ing of an optical soliton between two coupled fibers 1 8 1 yields 
such an example. 

Let us discuss the relevance of the various possible tunnel- 
ing effects lO. We mean now by tunneling din effect allowing 
a physical system in a double well potential to go from one 
well to the other either by purely quantum effects or by a mix- 
ture of thermal and quantum fluctuations or by thermal fluctu- 
ations only. This defines an effective tunneling time. The first 
typical time is the quantum tunneling time, Tq, then the co- 



herence time Tc. The latter time is the maximum time during 
which the quantum state inside one of the wells remains co- 
herent, supposing that the condensate is at finite temperature, 
since at zero temperature this coherence time is infinite. In the 
case of BEC, Tc is the mean-free flight time of a particle in the 
ground state colliding with a thermal particle divided by the 
number of particles in the ground state, the collision process 
being Poissonian. If there are many particles in the ground 
state, this coherence time may become very short because of 
this division. The third typical time is the Arrhenius-time Ta 
for the passage over the energy barrier separating the two clas- 
sical ground states under the effect of the thermal fluctuations. 
If Ta is very large, two possibilities are met. If Tq <^Tc, the 
effective tunneling time is Tq. If Tc <^ Tq, the effective 

tunneling time is of order ^ |[6l|. If the Arrhenius time is 
less than Tq, the transition will always be by activated barrier 

crossing. Otherwise, the shortest of the two times, ^ or Ta, 
will be the effective transition time. 

In the present problem the tunneling time can be calcu- 
lated by quantizing the dynamical GP equation |7|. If one 
does not attempt to compute the prefactor of the exponential 
ex.-p{SE/K), the derivation of Tq, that amounts to calculate 
the Euclidean trajectory, is usually fairly standard. However, 
because the present problem has a peculiarity compared to 
the ones usually treated, we shall explain the general idea in 
some details. The equation ^ is the Euler-Lagrange condi- 
tion making stationary the action 




drdtC, 



associated to the Lagrangian £, which is the following func- 
tional of ^ and ^, 

^ i f_^d^ -rd^\ 1,^.,2 1 ,.,4 

+ ^(F(r) - + ^ - . (4) 

We include the ground state energy /i in (|4]) to make it 
clearer the possibility of real or complex solutions. To per- 
form the analytic continuation of the dynamical equations, 
one has to write first the relevant quantities with real func- 
tions. Let us write ^as^ = ^/ + iv with u and v being real 
functions of r and t. The Lagrangian ^ becomes: 

where 

-Coo = ^ ((Vu)2 + (Vvf) + [vMiu) - uM{v)] 

+ ^{V{r)-^^){u' + v') + ^{u^+v^f. (5) 

In Feynman's formulation of quantum mechanics |5 1 any 
observable is given by a formal integral over all paths with 
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an integrand that is a function times an exponential ^ , 
S{X{-)) being the classical action of the path X{-). Here 
S{X{-)) would be calculated along all paths where the two 
functions u{t^r) and v{t^r) represent the vortex on one side 
(at time minus infinity) and the other with the vortex on the 
other side (at time plus infinity). Along the path, the wave- 
function ^(r, t)) = ix(r, t)-\-iv{T^ t) is not necessarily a solu- 
tion of the dynamical equation ([T]), but must tend at t = ±00 
to steady solutions of the GP equation (note that /i^ is sub- 
tracted in equation ([T]) to have a time independent solution). 
The tunneling corresponds to situations where the end-points 
of the trajectory X(-) are not linked by a classical trajectory, 
yielding the dominant saddle-point contribution to Feynman's 
integral at h small. In this limit, and if there is no such clas- 
sical trajectory, one extends the range of possible values of t 
to the complex plane. The Euclidean action associated to tun- 
neling is found by taking t purely imaginary in the equations 
of motion, t = ir,r real, leading to the expression 

= i (^^^ ~ ^^7) ^ 

The Euclidean equations of motion for u and v are derived 
by variation of J dt J (IyCe, leading to the two coupled equa- 
tions: 

= -\^^^ + + + + ^ - /i), (7) 

and 

i^ = \v^u-M{v)-\{u' + v^ + V (8) 
or 2 

If the contribution proportional to M is absent (that is with- 
out rotation), the Lagrangian in equation ([5]), and the dynami- 
cal equations (|7])-([8]) become real by changing v into iv, yield- 
ing a formally real Euclidean dynamics. However this can- 
not be done with a non-zero M because the transformed La- 
grangian is complex and cannot be made real by any simple 
transformation. Although this is rarely considered, a fully 
complex Euclidean action does not hurt any general principle. 

In the present problem, the tunneling probability is propor- 
tional to exp(3?(S'£;)/ h), and the end states, at times r = ±00, 
are the classical ground states, with real u^o and ^oo, because 
either in real or in Euclidean dynamics, the equilibrium state is 
the state which minimizes J drCoo (i.e. the part of the action 
not proportional to time derivatives). Note that by substitut- 
ing t for —ir in the pair of equations ^ and ([5]), one recovers 
the classical GP equation, where u and v are purely real. In 
the case considered here, where the asymmetrical one- vortex 
state is the ground state, the GP equation leads to stable equi- 
librium states, (Uoo and Voo), that forbids any tunneling. While 
the tunneling occurs when using the Euclidean equations 
(jSj with real r, where (Uoo and Voo) are unstable equilibrium 



states. For arbitrary r, the functions u and v are both com- 
plex, without any simple relation in between, like being a pair 
of complex conjugates. Therefore the present Euclidean prob- 
lem involves two complex fields, each having their own set of 
different vortices, a rather unusual fact. The numerical inte- 
gration of the Eulerian system requires some caution, and will 
be treated in a future publication. Actually, while the math- 
ematical problem is well posed, as a Dirichlet problem with 
boundary conditions at infinite time, the numerical search of 
solutions should avoid to solve an initial value Cauchy prob- 
lem that is ill-defined because of the instability of this Eu- 
clidean dynamics at very short wavelength |9|. Furthermore, 
a change of variable r into, for example tanh(r), would lead 
to a finite integration time. Then the trajectory joining the two 
equilibrium states, could be found by starting slightly off one 
equilibrium state and by finding the pair v that minimizes 
I J drdrjC{u, v)\'^. This procedure should avoid the instabil- 
ity effects for large wave numbers. Along this trajectory, the 
transient vortices formed on each field, will merge at the end, 
into a single one. 
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